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\ The arn^/my^ order supersymmetric electroweak corrections arising from 

loops of chargino, neutralino, and squark to top quark pair production by 
gg fusion at LHC are calculated in the minimal supersymmetric model. We 
Q_i' found that the corrections amount about a few percent. 
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I. Introduction 

The top quark has been found experimentally by the CDF and DO Collaborations with 
the mass and production cross section nit = 176 ± 8{stat) ± 10{syst) GeV a = GStl'lpb, 
and mt = 199lLll{stat) ± 22{syst) GeV (x = 6.4 ± 2.2p6, respectively 0. Although this 
measured mass is close to the central value predicted by the best fit of the Standard Model 
(SM ) to the latest LEP data, the central value of the cross section is somewhat larger than 
the Standard Model prediction ati = 5.52lo;45 pb for rrit = 175 GeV at ^/s = 1.8 TeV pp 
collider in which the effects of multiple soft-gluon emissions have been properly resummed 0. 
In addition, there are still a number of unsolved theoretical problems in the SM. New physics 
beyond the SM are still possible. Among various models of new physics so far considered, 
supersymmetry (SUSY) is a promising one at present. The simplest and interesting SUSY 
model is the minimal supersymmetric extension of the standard model (MSSM) At 
the future multi-TeV proton colliders such as the CERN Large Hadron Collider (LHC), ti 
production will be enormously larger than the Tevatron rates and the accuracy with which 
the top quark production cross section can be measured will be much better (the uncertainty 
is about 5% at LHC [Q). Thus theoretical calculations of the radiative corrections to the 
production of the top quark at those colliders are of importance. 

QCD corrections to 0{a^) and electroweak one- loop corrections in the SM to tt produc- 
tion in hadron colliders are carried out in Ref. [0 ||5[ and Ref. [^. Yukawa corrections to ti 
production at the Fermilab Tevatron and LHC in two-Higgs-doublet models are calculated 
in p. In the MSSM, electroweak corrections from chargino, neutralino and squark to the 
top pair production via qq annihilation at the Fermilab Tevatron are calculated in Ref. 0. 
Recent calculation of the supersymmetric QCD corrections to the top quark pair production 
at the Tevatron shows that they increase the cross section by about 20% [[l^] [|lT|. But this 
is still within the experimental uncertainty 30%. At LHC the main production mechanism 
of top quark pair is the gluon-gluon fusion process gg ti. In this paper we investigate the 
electroweak corrections of order amf/m'^ arising from chargino, neutralino and squark to 
the top quark production by the process gg — > ti at LHC. The formalism of the calculation 
of the corrections to the matrix elements will be given in Sec. II. In Sec HI, we present our 
numerical examples and discussions of the corrections to the cross sections in the MSSM. 

II. Formalism 

The tree-level Feynman diagrams and the relevant supersymmetric electroweak correc- 
tions to gg — >■ ti are shown in Fig.l (u-channel diagrams of (b) and (e)-(h) are not explicitly 
shown) in which the dashed lines in the loop represent the squark t^ or bi {i = 1, 2) and the 
solid lines represent neutralinos or charginos, respectively. 

The supersymmetric partner of left- and right-handed massive quarks mix [|l^. The 
mass eigenstates qi and q2 are related to the current eigenstates and 

qi = qLCOs9g + qRsmOg, q2 = -qismOg + qRCOsOg (1) 

The mixing angle 6t and the masses m^^, m^^ can be calculated by diagonalizing the 
following mass matrix 
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* \mtmLR j 

Ml = m?^ + + - ^ sin^ Qw) cos(2/?)m| 
2 

= m?^ +m] + - sin^ 6*11/ cos(2/5)m| 
"^Li? = cot /5 + A (2) 

where m?^ , m?^ are the soft SUSY-breaking mass terms of left- and right-handed stops, /x is 
the coefficient of the Hi H2 mixing term in the superpotential, At is the parameter describing 
the strength of soft SUSY-breaking trihnear scalar interaction tLtRH2, tan/5 = V2/V1 is the 
ratio of the vacuum expectation values of the two Higgs doublets. 
From Eq.(|D, we can get the expressions for and 6t : 



*1,2 



2 + M4 t ^{Ml - Mry + 4m2my (3) 
i{met = ^ ^ (4) 

For the sbottoms, we neglect the mixing between the left- and right- handed 
sbottoms(6'b = 0) and have 

^li,2 = ^h,bR + "^6 ± {TI r - Qb sin ew) cos(2/3)m|, (5) 

where T|^^ = 0, Qb = — | and m?^ are the soft SUSY-breaking mass terms for left- 
and right-handed sbottoms. 

In the presence of squark mixing,the squark-quark-neutralino and squark- quark- chargino 
interaction Lagrangian of order grrit/mw is given by 
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where g is the SU(2) coupling constant , and aj^j, b^^j, aj^j, bf^j, ai^-, b^^., ai^-, bi^. are 
given by 

'^kj = ^ (^i4 cos Ot + Nj4 sin Ot), b^^j = ^ {N*^^ cos 6t - iVj4 sin Ot), 

ai^j = ^ (-^^4 sin Ot + Nj^ cos 9t), b^^j = ^ (-^4 sin 9t - Nj^ cos 9t), 

"feii = Kj = cos ^fe, = bi^^ = -Wj^ sin Oh (7) 

Vj2 are the elements of 2 x 2 matrix V and Nj4^ are the elements of 4 x 4 matrix N (see the 
Appendix). 
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At the tree level, the S-matrix element is composed of three different production 
channels(s-,t-,u-channel) as follows: 

= -liT^'T^ - T^T'')jiM^', (8) 
M* = -zg',{T'T^),MP2)Mi + m)hv{pi)/{i - ml) 

= -z{T'T^),,M',', (9) 

= -i{T-T').j,Mi;\ (10) 

where q = P2 — Pi, £4 = ^'^{Pi), £3 = ^^{Ps), is given in the Appendix. Instead 
of calculating the square of the amplitudes explicitly, we calculate the helicity amplitudes 
numerically by using the method of Ref . [0 . This method greatly simplifies our calculations. 

The 0{aml/myy) SUSY electroweak corrections to gg tt are shown in Fig.l (c)- 
Fig.l (h). The sum of them is QCD gauge invariant without the strong coupling constant 
renormalization . In our calculation, we use dimensional regularization to regulate the 
ultraviolet divergences and adopt the on-mass-shell renormalization scheme. We also discard 
the terms proportional to 75. 

We only give the explicitly results of the s- and t-channel contributions to the SUSY 
electroweak corrections. The u-channel results can be obtained by the following substitu- 
tions: 

P3^P4, T'^^T^ i^u. (11) 

Fig. 1(c) lead to the s-channel vertex correction ^M**^: 

6M^' = -tgl{zU,T%u{p,){F^' + F^'f + ^Xpi)/^ 

= -i{T^T^ - T^T'')5M'^' (12) 

Fig. 1(d) gives 5M'^: 

6M^' = -zgliT^T' + T''T%MP2)Ffv{p,) 

= -{{T^T" + T^T^)ji5M''^' (13) 

The top quark self-energy 5M'^^''^'^ of Fig. 1(e) is: 

= -i{T^T^)ji6M'^^^'" (14) 

Vertex correction 5M^^'^ of Fig. 1(f) is: 

5M^i'* = -^g',iT'^T%MP2)iF;;''' + Ff^^ + rt'M + m,)i,v{p{) / {t - m?) 

= -i{T^T'')ji6M''^''' (15) 



v2,t 



Vertex correction 6M of Fig. 1(g) is: 
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= -i{T^T'')ji6M''^'" (16) 
^^box,t ^Yie box diagram Fig. 1(h) is: 

5M^-'* = -^gliT'T'^),^u{p,)[F'/ + /^a^'Xpi) 

= -i{T''T^)ji6M^'"''" (17) 

The total amphtude can be written as: 



M,, = -z[{M+ + 6M+)0^P + (M, + (5M-)0], ^ ], (18) 



where 



M(+ = M*' + Mo"' , Mq = M*' - Mg"' - 2M^', 

6M- = -26M'^' + 6M'^^^'" - 5M^^'^'"' + 5M^^'*' - 5M^^'"' 

+(5M^^'*' - 5M^2'"' + 5M^°^'*' - 5M''°^'"' (19) 

The color sum of the corrected amplitude square is: 
^ |M|2 = I|Mo+p + 3|Mo-|2 

color 

+ ^i?e(Mo+5M+t) + 6i?e(Mo"5M~^) (20) 

The spin sum as well as phase-space integration and parton distribution convolution are 
done by the VEGAS program. The correction cross section Aa is defined as 

Aa = a - ctq, (21) 

where a is the cross section given by |Mp and ctq is the tree level QCD cross section given 
by |Mo+P and \Mq\'^ . 

III. Numerical Examples and Discussions 

The production cross section is obtained by convoluting the partonic cross section with 
certain parton distribution. The relative correction is not sensitive to parton distribution. 



In this paper, we take the Martin-Roberts-Stirling (MRS) parton distribution set A' |jT4 
with = s. The following parameters are used in our calculation: 

= 14 TeV, mt = 176 GeV, rrib = 4.9 GeV, 

mw = 80.22 GeV, mz = 91.175 GeV, a = — (22) 

128.8 
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Care must be taken in the calculation of the form factors expressed in terms of the standard 
loop integrals defined in Ref. |15|. As has been discussed in Ref. the formulae for the 



form factors given in terms of the tensor loop integrals will be ill-defined when the scattering 
is forwards or backwards wherein the Gram determinants of some matrices vanish and thus 
their inverses do not exist. This problem can be avoided by taking the kinematic cuts on 
the rapidity y and the transverse momentum px- In this paper, we take 

\y\ < 2.5, pt > 20 GeV. (23) 

The cuts will also increase the relative corrections 0. 

We first checked the QCD gauge invariance by the substitution p4 — > 64 and pa £3 and 
find that 6M~^, 6M~ are a few order of magnitudes smaller. 

In the calculation of the chargino and neutralino masses, we fix M = 200 GeV, fi = 
— 100 GeV and use the relation M' = ^2M{see the Appendix). We also assume m^-^ = 

The relative correction to the hadronic cross section as a function of the squark mass 
parameter rriq with tan j3 = 1 and rriLji = (corresponding to non-mixing case) is presented in 
Fig. 2. For tan/5 = 1, the chargino masses m-+ = (220, 120) GeV and the neutralino masses 

m^o = (105,221,128,100) GeV. The correction is always negative. For rriq < 150 GeV, 
the correction is very sensitive to rrig. A sharp dip at about rrig = 56 GeV is due to the 
singularity of the top quark wave function renormalization constant at the threshold point 
rrit = + m-+ (note that for tan/5 = 1, m^,^ = nig). This singularity will disappear if 
the finite widths of the top quark and the charginos are taken into account. The correction 
exceeds -5% only in a small region near the dip. The correction approaches to zero at large 
rriq which shows the decoupling behaviour. 

Fig. 3 shows the dependence of the relative correction to the hadronic cross section on the 
stop mixing parameter rriLR. We set rriq = 100 GeV and tan (3 = 1. rriLR affects the mass 
splitting and mixing angle of ti and ^2- The mass splitting increases as miR increases. We fix 
ti to be the light one(cf. Eq.(3)). Therefore, 6t = \ for miji < 0, 9t = — | for m^ji > 0. The 
mixing angle causes the asymmetry of the relative correction between > and rriLji < 
although the mass splitting of ii and ^2 is symmetry between rriLR > and miji < 0. The 
dip at about miji = —200 GeV is due to the threshold effect at rrit = mi + m^o. No dip 

is found at rriLR = 200 GeV because when rriLR > 0, 9t = — 7r/4, the ttix'j coupling is 
proportional to 75 {aij = 0, bij 7^ 0). From the expression of the top quark renormalization 
constant (see the Appendix), one can see that the singularities of Go and Gi cancel with 
each other when a^j = 0, bij 7^ 0. 

In Fig. 4, we present the tan (3 dependence of the relative correction to the hadronic cross 
section at given rriq = 100 GeV and miFt = 100 GeV. tan/? slightly affects the stop mass 
splitting and . The factor 1/ sin^/? in the coupling constant leads to the rapid increase 
of the correction in the range tan/5 < 1. But the increase is somewhat more quickly than 
1/ sin^/3 because m^^ decreases as /5 decreases. 

From Fig.2-4, we see that only for tan/? < 1 and a small region near the threshold t — 
bix^ and t — > tix^ the correction may exceed —5%. Otherwise, the correction amounts only a 
few percent smaller than —5%. Therefore, we conclude that the supersymmetric electroweak 
corrections of order amf/m^r to top quark pair production at LHC are potentially observable 
for tan/5 < 1 and small parameter region near the threshold t —>■ biXj , t tix^ . 
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Appendix 

We give here the form factors for the matrix element appeared in the text. They are 
written in terms of the conventional one-, two-, three- and four-point scalar loop integrals 
defined in Ref . . 

^o' = - E E _ cPL[mA^l - bl)m - pi) ■ rco 

3 i=ii,i2,bi,b2 

-2C^° ■ r)](-p2, k,mj,mi,mi) 



Fi' = E _ E _ CPL[6Z^^] 

3 i=ii,i2,bi,b2 

Ft' = - E E _ CPLiial + hl)m - pi) . VC^'^ 

3 i=ii ,i2 ,bi ,b2 



= E CPL[mt{al + bl)Cn - mj{al - bl)Co]{-p2, k, m^, m^, m^) 

3 i=ii,i2,bi,b2 

3 j=?l,?2 ,61,62 
3 i=ii ,i2 M ,b2 

Fq^'^ = - E E CPL[2mj{al - hjj){p2 ■ e^Co - ■ e4)](-p2,P4, rrij, rrii, mi) 

3 i=ii ,i2 ,bi ,b2 

Ft' = E E _ CPL[6Z^,] 

3 i=ti ,t2 ,bi ,b2 



Ft''' = -E E _ CPmal + bl){p2-e,C'''^ 

3 i=ii ,i2 ,bi ,b2 

-C^^''(e4))](-p2, P4, rrij, m,, m,) 



Ft' = - E _ E _ CPL[2m,{al - l^,){-p, ■ e^Co 

3 i=ii M ,bi ,b2 

• e3)]{pi, -p3,mj,mi,mi) 
Ff' = Y: _E _ CPL[5Zl] 

3 i=ti,t2,b\,b2 

Ff'' = -Y: E _ CPL[2{al + hl){-p,-e,C''>^ 

3 i=ti ,t2 ,bi ,b2 

-C^^''(e3))](pi, -p3,mj,mi,mi) 



Ft = E E CPL[Am,{al - bl){p2 ■ ■ e^Do - Pi ■ e,D^^ ■ 64 

3 i=ti,t2,b\,b2 

+P2 ■ 64^^° ■ 63 - £3 ■ D'^^{e^))]{-p2,P4,P3,mj,mi,mi,mi) 
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i=ti,t2,bi,ti2 

+P2 ■ €4-0^^^(63) - D^^>'{e3,e4))]{-p2,P4,P3,mj,mi,mi,mi) 

In the above i, j summation, j = 1, 2, 3, 4, rrij = m^o for i = ti, t2 , j = 1,2, rrij = 
for i = 61, 62- As 6*6 = 0, 62 actually does not contribute to the sum(cf. Eq.(|^). 

k = pi+p2=pz+PA, s = k'^, i = q"^ = {p2 - Pif, U = {p2-P3Y, 
r'^ = (-P4 + PsTe^ ■ £4 + (2p4 + Ps) ■ 6364 - (2P3 + P4) ■ e^e^ , 



CPL 



am 



Svrm^ sin 9w sin /? 



In our calculation, we calculate the tensor loop integrals C^, C^'^ , D'^, and D'^""' 
numerically instead of expanding them explicitly. 
The renormalization constants are: 



5m. 



[{oHj + ^%){Bi + 2mlGi) + 2(a.j. - bl)mtmjGo]{p^, m], m- )|p2=^2 



dBo{p'^,mi,mi) _ dBi{p^ ,m[,mi 

— , (jri — — 



where Gq ■ , „ , 

The neutralino and chargino masses , the 2x2 matrix V and the 4x4 matrix are 
given by the following relations 0: 



X 



M mwV2si-n.(3 
mw\/2cos (3 ji 



D 



Y 



M' 






M 



—mzsw cos P mzsw sin P \ 
mzCw cos P —mzCw sin P 
mzSw cos P mzCw cos P — /i 

V mzSw sin /? — m^cvi/ sin P —ji 



D 



where M, M' are the masses of gauginos corresponding to SU{2) and U{1), respectively. 
With the grand unification assumption, we have the relation M' = ^[g''^/g'^)M. M|, is a 
diagonal matrix and the chargino mass squares correspond to its diagonal elements. 
is a diagonal matrix with non-negative entries which give the masses of neutralinos. The 
chargino and neutralino masses depend on the parameters M, yU and tan p. 
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Figure Captions 

Fig.l Feynman diagrams of tree level and 0{am^/m^r) SUSY electroweak correction. 
Fig. 2 Relative correction to hadronic cross section versus rriq with tan P — 1 and rriLR — 0. 
Fig. 3 Relative correction to hadronic cross section versus with tan (3 — 1 and rrig — 100 
GeV. 

Fig. 4 Relative correction to hadronic cross section versus tan /? with rriq — 100 GeV and 
rriLR = 100 GeV. 
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